/9/n & s 

NASA CONTRACTOR REPORT 177487 

/S7.&9J 4 - 

p-US- 


VSAERO Analysis of Tip 
for the Tree-Tip Rotor 


.Planforms 


D. M. Martin 
P- E. Fortin 


(EASA-CB-177487) VSAERO ANALYSIS OF TIP 
PIANFORHS FOB TEE FBEI-TIP EOICB (Kansas 
Uciv. Center for Research) 125 p CSCX 01A 


V 


N88-28036 


Unclas 

G3/02 01570S9 


CONTRACT ncc 2-175 
June 1988 


m/vsA 



NASA CONTRACTOR REPORT 177487 


VSAERO Analysis of Tip Planforms 
for the Free-Tip Rotor 


D. M. Martin 
P. E. Fortin 

The University of Kansas Center 
for Research, Inc. 

Lawrence, Kansas 


Prepared for 
Ames Research Center 
under Contract NCC 2-175 
June 1988 


NASA 

National Aeronautics and 
Space Administration 

Ames Research Center 

Moffett Field, California 94035 



ABSTRACT 


This report presents the results of a numerical analysis of two 
interacting lifting surfaces separated in the spanwise direction by 
a narrow gap. The configuration consists of a semispan wing with 
the last 32% of the span structurally separated from the inboard 
section. The angle of attack of the outboard section is set 
independently from that of the inboard section. In the present 
study, the three-dimensional panel code VSAERO is used to perform 
the analysis. Computed values of tip surface lift and pitching 
moment coefficients are correlated with experimental data to 
determine the proper appraoch to model the gap region between the 
surfaces. Pitching moment data for various tip planforms are also 
presented to show how the variation of tip pitching moment with 
angle of attack may be increased easily in incompressible flow. 
Calculated three-dimensional characteristics in compressible flow at 
Mach numbers of 0.5 and 0.7 are presented for new tip planform 
designs. An analysis of sectional aerodynamic center shift as a 
function of Mach number is also included for a representative tip 
planform. It is also shown that the induced drag of the tip surface 
is reduced for negative incidence angles relative to the inboard 
section. The results indicate that this local drag reduction 
overcomes the associated increase in wing induced drag at high wing 
lift coefficients. 
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1. INTRO DUCTION 

I • i FREE— TIP CONCEPT PROJECT STATUS 

The free-tip rotor concept has been the subject of a 
considerable amount of research in the past three years. Its main 
purpose is to improve rotor performance and reduce vibration. The 
free-tip is structurally decoupled from the inboard portion of the 
blade and allowed to pitch about an axis located at the quarter 
chord of the inboard blade (Figure 1.1); a positive control moment, 
opposing the negative moment due to lift, is applied to the tip via 
the pitch shaft. As the blade rotates around the azimuth, the tip 
weathervanes about an equilibrium position, thus yielding finite 
lift. 

Since helicopter forward flight is characterized by large 
radial and azimuthal variations in blade loading, the relatively 
uniform lift produced by the free-tip should improve the overall 
blade lift-over-drag ratio (L/D). This is achieved by eliminating 
the negative tip lift on the advancing side usually associated with 
flight at high advance ratios. Since the amplitude of the lift 
oscillations should also be reduced, the free-tip equipped rotor 
will have improved vibration characteristics. 

To date, the aerodynamic characteristics of free-tip 
configurations have been studied experimentally and analytically. 

In early 1982, an investigation of tip planform influence on the 
aerodynamic loading of an unswept wing was conducted in the 
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Figure 1.1: Free-tip configuration 






NASA 7xl0-foot wind tunnel (Reference 1). Several tip 
configurations were tested and analyzed; this led to the selection 
of a tip geometry which yielded the desired aerodynamic 
characteristics. The next phase in the experimental program, which 
took place in November 1983 (Reference 2), consisted of the testing 
of a small-scale free-tip rotor. Each of the four blades was fitted 
with a free-tip which was attached to the inboard portion at the 1/8 
chord point (Figure 1.2). The tip planform geometry was based on 
the results of the 1982 semi-span test. The passive controller 
generated a nose-up pitching moment which varied with rotor RPM. In 
the course of this test, the free-tip 's response was measured for 
advance ratios from 0.1 to 0.397, and some performance gains were 
observed. Another semi-span test was conducted in January 1985: 
the spanwise variation of wing loading was measured by computing the 
circulation at various spanwise stations from Laser Doppler 
Velocimeter data. Flow visualization studies were also performed to 
determine the geometry and behavior of the strong vortex which is 
shed at the tip attachment point. A formal report on this latest 
wind tunnel investigation is still pending. The ongoing analytical 
effort mainly consists of analyzing the free-tip concept using 
advanced computer codes. 

1.2 FREE-TIP RESPONSE 

At this stage in the experimental and analytical study of free- 
tip performance, the aerodynamic criterion used to measure tip 
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Figure 1*2: Free-tip-equipped rotor in wind tunnel 


response Is the pitching moment derivative (C m )j. The magnitude of 

a 

this parameter is strongly dependent on the offset between the 
aerodynamic center of the tip planform and the pitch axis. If the 
tip is swept back, the aerodynamic center is moved aft; with the 
lift acting at the tip aerodynamic center, the total aerodynamic 
moment (excluding the controller moment and drag force) about the 
pitch axis may be written (Figure 1.3) 

M A = k a a = [(C L >T a< l S tip jAc + c m qS tip' (ia > 

a r ac 

where k a is the resultant aerodynamic spring rate and Ac is the 
offset between the pitch axis and the aerodynamic center and is 
usually negative. Note that a here is measured from the zero-lift 
line. The aerodynamic spring. rate, in pounds-feet per degree of 
angle of attack, may thus be expressed as 

k a - 1(C L VW* + C m <‘- 2 > 

a r ac r 

The pitching moment derivative (C ) T is usually expressed in terms 

m a 1 

of an arbitrarily selected reference point (which for the free-tip 
is the shaft pitch axis), 

<«. >T ■ < C L XAC ) * <°L < 1 ' 3 > 

a a c a 

then Equation (1.2) may be rewritten: 

k « - I(C » >T + C „ '“HW <l ' 4 > 

a ac 

In Reference 3 It is shown that k Q appears explicitly in the 
equations of motion of the free-tip. As one of the parameters in 
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the coefficient of the pitch degree of freedom (first order term), 

k Q also appears in the expression for undamped natural frequency. 

Clearly, from Equation (1.4) and its relationship to natural 

frequency, the derivative (C m ) T has a significant impact on the 

m ct 

response of the free-tlp in steady flow and is thus used as the 
primary design parameter in this study. 


1.3 USE OF THE PANEL CODE VSAERO 


Because of the technical difficulties associated with the 
acquisition of tip lift and pitching moment data from an actual 
rotor in a wind tunnel, much of the aerodynamic data accumulated so 
far were obtained from tests of fixed semi-span wings with 
structurally decoupled tips. The tips were not allowed to pitch 
freely but were set at specified incidence angles with respect to 
the inboard section. The aerodynamic characteristics of the tip 
(C; , C , and C n ) were then obtained from force and moment data 

Lh£* THip 1 U<p 

measured by a balance at the tip junction. 

Obviously, to obtain the best possible prediction of free-tip 
aerodynamic parameters, it is necessary to account for the Induced 
effects of the Inboard wlng/blade on the tip. Since the loading of 
a helicopter blade is generally higher at the tip (because of the 
radial variation of dynamic pressure), the semi-span configuration, 
with its near half-elliptical spanwise loading distribution, does 
not adequately represent the physical problem. Therefore, the 
crucial assumption which has been made in the past and continues to 
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apply in the present study is that there exists a relationship 
between the aerodynamic characteristics of a decoupled wing tip in a 
semi-span configuration and those of a free tip installed on a 
rotating helicopter blade in forward flight* 

As a result of this assumption, analytical tools which were 
originally developed for the study of fixed-wing configurations may 
be used to predict the parameters of interest for free-tip 
applications. The computer program VSAERO, a code designed for 
calculating the subsonic aerodynamic characteristics of arbitrary 
configurations, lends itself well to this task. The program allows 
for the modeling of separate lifting surfaces within a configuration 
and accounts for the aerodynamic interaction between them. Wake 
relaxation can also be accomplished by automatic updating of the 
orientation of wake panels through an iterative process (see Figure 
1.4). The potential flow solution is also coupled with a boundary 
layer routine; in subsequent viscous/potential iterations, the 
potential flow boundary conditions are modified by a source 
transpiration technique. Corrections for subcritical 
compressibility effects are also possible; as will be discussed 
later, this feature was used in the present analysis. 

Although the code provides a method to define regions of 
separated flow for the computation of nonlinear effects, this 
capability was not used. For the range of incidence angles in which 
the free tip typically operates, the flow separates on the lower 
surface. Several attempts were made to "stitch" the wake along a 
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Figure 1.4: VSAERO program global flowchart 












separation line defined from boundary layer data, but the computed 
aerodynamic parameters were highly unrealistic. This feature of 
VSAERO is still in the development phase at Analytical Methods, 

Inc., the originator of the program. 

1.4 OBJECTIVES OF THE CURRENT STUDY 

The objectives of the current study were twofold. First, the 
VSAERO program was validated for this application by comparing the 
aerodynamic parameters obtained from various modeling methods with 
the experimental data of Reference 1. This code was selected 
because it offered a wide range of possibilities in panel density 
selection, component identification and modeling, and number of wake 
shape and viscous/potential iterations. 

Second, a new family of tip planforms was designed and 

analyzed. This was required to increase the magnitude of (C ) T of 

m a 

the baseline configuration, and thus achieve better response of the 
tip to angle-of -attack variations. A steeper pitching moment slope 
is obtained by increasing the sweep angle of the leading edges: 
this moves the aerodynamic center aft, increasing the moment arm to 
the pitch axis. Although a target range of -0.012/deg for (C m 
was originally selected, intermediate geometries are also required, 
since most of the terms which play an important role in the free-tip 
equations of motion have not been set for a final design yet. These 
include, for example, the free-tip inertia about the pitch axis and 
the restoring moment applied by the controller. The aerodynamic 
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characteristics of the new planforms will serve as an initial data 
base for future free-tip designs. The lift and pitching moment 
behavior of the new tips and the semi-span wing drag characteristics 
due to each configuration were also investigated. 

The basic geometry used in this analysis is that of the fixed- 
wing semi-span configuration described in Reference 1. 
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2m CORRELATION STUDIES USING VSAERO 


This section describes the approach used to model the semi-span 
wing and tip using VSAERO. The predicted lift, pitching moment and 
drag of two tip planforms are compared to experimental data; these 
correlations are analyzed in detail to determine the best modeling 
method for the design and analysis of new tip planforms. 

The theoretical formulation of the VSAERO code is beyond the 
scope of this report. However, a summary of the analytical 
development of the program may be found in Appendix A. A more 
comprehensive description of the solution procedure is given in 
Reference 4. 

2 . 1 PRELIMINARY CONSIDERATIONS 

Before proceeding with the analysis of any configuration with a 
program of this nature, it is necessary to select the optimum number 
of body and wake panels along with the required number of wake shape 
and viscous/potential iterations. The configuration which was 
considered to determine these parameters consisted of two distinct, 
separate surfaces with the tip moderately swept (Reference 1); in 
previous studies, this particular modeling approach was preferred 
because it closely resembled the actual configuration. 

The method which was applied to determine the best panel 
arrangement is relatively simple. First, a full-cosine distribution 
was selected for chordwise panel spacing: this results in a higher 
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panel density at the trailing and leading edges and is desirable for 
good chordwise loading prediction. A half-cosine distribution 
defining the spanwise regions was selected for the inboard blade and 
the inboard and outboard halves of the free tip. This approach 
generates more panels in the vicinity of the structural 
discontinuity and the outboard extremity of the free tip which are 
regions of large spanwise loading gradients for certain tip 
incidences. Once these parameters were set, the smallest number of 
panels was determined by increasing their number progressively until 
convergence of C L and C m values for both lifting surfaces was 
achieved. 

The final combination therefore consists of 15 chordwise 
divisions (on both the upper and lower surfaces) for all airfoil 
sections defined in the configuration; the inboard blade is divided 
into 7 spanwise columns while the tip is divided into 6 similar 
regions. The body surface therefore consists of 495 panels, as 
shown in Figure 2.1. This number Includes 15 vertical panels 
closing each of the tips of both components at the gap. As will be 
seen later, this approach was reevaluated in the correlation 
studies. There are 75 vertical panels on the most outboard tip 
section. 

The spanwise wake panel distribution is automatically set by 
the selected number of wing panel columns; i.e., each column of wing 
panels sheds one column of wake panels. The streamwise number of 
divisions in the wake was set to 20, with a higher panel density at 
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Figure 2.1 (a): Semi-span wing configuration with TIP 1A 
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Figure 2.1 (b): Semi-span wing configuration with TIP 1A (top view) 



the wing trailing edge. The total number of wake panels therefore 
amounts to 260 (13 columns x 20 divisions, Figure 2.2). However, 
since the upstream edge of the wake was defined (by the present 
user) slightly behind the quarter chord line, some wake panels 
overlap the body and are not shown. The upstream edge was selected 
at this location because it was originally intended to stitch the 
flow separation line. Unfortunately, the use of this option was not 
successful, due to lack of information and limited user 
experience. Therefore, for all cases compiled in this report, 
separation was assumed to take place at the trailing edge. Even if 
some wake panels are tentatively defined upstream of the trailing 
edge by the user, the code will correct automatically if trailing- 
edge separation is selected later in the wake geometry definition. 
For this reason, only 182 wake panels are actually needed, as shown 
in Figure 2.2. The downstream edge of the wake is located roughly 
7.5 chord lengths behind the trailing edge. 

Generally, the solution converged for a small number of 
viscous/potential iterations. Since only the right-hand sides of 
the system of linear- equations to be solved by the code (see 
Appendix A, Equation A. 14) are affected by the source transpiration 
values returned from the boundary layer routine, specifying 
viscous/potential iterations does not result in prohibitive CPU 
costs. In the present application, three such iterations were 
considered to be adequate. When the wake relaxation loop Is 
executed, the wake panels are repositioned so that their streamwise 
edges are aligned with the local calculated flow direction. 
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Figure 2.2 (b): Semi-span wing and tip wake panels (top view) 



Equations (A. 14) also illustrate the fact that wake shape iterations 
require that the wake influence coefficients be recomputed, 
increasing the required CPU time. Therefore, one wake shape 
iteration per viscous/potential iteration was used as 'a compromise 
to reduce job processing times. 

Two integral boundary- layer routines are available in the 
VSAERO program. The first is based on the computation of boundary- 
layer characteristics along external surface streamlines. The 
procedure accounts for surface curvature and streamline 
convergence/divergence under the assumption of local axisymmetric 
flow. In earlier studies, when the leading edges of the free-tip 
were only moderately swept, the procedure was acceptable. However, 
Reference 4 shows that this 2-D method could break down in regions 
of high crossflow. Recent investigations with swept tips have 
confirmed this finding. In cases where some regions of separated 
flow (however small) were encountered, the computed local skin 
friction drag contribution was large and negative. In view of this 
problem, the second boundary-layer routine, which includes a 
crossflow model, was used for the analysis of the new 
configurations. The user prescribes a series of chordwise strips 
across the wing surfaces where the analysis is desired. 

The preceding figures (2.1 and 2.2) illustrate the final body 

\ 

and wake panel distributions which were used for the correlation of 
tip lift, drag, and pitching moment. The wind tunnel walls and test 
section were not modeled for the correlation phase of this study. 
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2.2 VS AERO CORRELATION WITH WIND TUNNEL EXPERIMENT 


The main purpose of this phase of the research is to obtain the 
best possible correlation of and with experimental 
results. Also, an investigation of the tip’s drag polars will be 
carried out. First, it is necessary to define the tip incidence 
relative to the inboard blade as 


A9 = (2.1) 

where the angle of attack in the tests of Reference 1 was measured 
relative to the tunnel floor. In this numerical application, 
and are measured from the x axis of the global coordinate 
system (which is located at the quarter chord of the section at the 
plane of symmetry); the free-stream velocity vector is parallel to 
this axis. From the rotor test results of Reference 2, the A9 of 
the tip is usually negative for all azimuthal positions of the 
blade: however, in the semi-span wing experimental results of 
Reference 1, the A9 f s were limited to values of +5.0, 0.0, and -5.0 
degrees. Since the data of Reference 2 represent more recent 
findings on free-tip behavior in forward flight, it was decided to 
limit the correlation to the only negative value of A9 for which 
experimental data are available, i.e. A0 » -5°. 

The airfoil section used for all free-tip designs discussed in 
this report is the V23010-1.58, a modified version of the NACA 23010 
section. Its most noticeable feature is the drooped leading edge 
region. The twist distribution of the inboard wing is nonlinear, 
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with a maximum variation of 0.321 degrees from the root to the 
beginning of the tip section. The free-tip planforms considered 
have all been designed without twist. More information on the 
V23010-1.58 airfoil and on the semi-span twist distribution is 
contained in Reference 1. For the correlations that follow, two 
configurations taken from Reference 1 were modeled using VSAERO. 
The planform parameters appear in the following table: 




Table 2.1: 

Geometric 

Parameters 

for TIP 

IA and TIP 

IB 



Ale 

A c / 4 

X 

C INB 

c 

s t 



(deg) 


— 


On) 

(m 2 ) 

TIP 

1A 

42 

35 

0.3 

0.208 

0.170 

0.0532 

TIP 

IB 

48 

45 

0.6 

0.208 

0.187 

0.0683 


Both tip planforms were analyzed for correlation of lift, drag, and 
pitching moment. In this chapter, the computed values of pitching 
moment coefficient are normalized by the mean geometric chord of the 
tip surface (c). 

2.2.1 Lift and Pitching Moment Correlation 

Three different methods were used to model the structural dis- 
continuity between the outboard surface and the Inboard section: 

• two separate components, closed tips 

• two separate components , open tips 

• single continuous surface, with a step change in incidence 
at the tip junction 
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Figures 2.3 and 2.4 illustrate the tip paneling arrangement for 
the first method, which initially seemed to be the most logical 
approach. The two surfaces are separated by a gap of roughly 3 mm; 
observe that both the inboard section tip and most inboard extremity 
of the free-tip have been paneled. This is the most distinctive 
feature of this particular approach. The results obtained from 
VSAERO are plotted against experimental data in Figure 2.5. 

To obtain these data, the angle of attack of the entire 
configuration (wing and tip) was varied with A0 =* -5° held constant 
as discussed earlier. The results indicate that the predicted tip 
lift is nonlinear for the entire range of under consideration. 
This effect is believed to be related to potential flow effects in 
the gap region. Since the boundary-layer method which accounts for 
crossflow effects was used in this study, viscous effects could not 
be modeled in the gap between the surfaces. Therefore, the fluid 
flows through the narrow gap and around the corners created by the 
end cap panels and upper and lower surface panels. This results in 
high local velocities and large values of Cp on the panels adjacent 
to the gap. Another possible factor in this problem could be that 
of control point/panel edge proximity. Since a doublet panel is 
equivalent to a quadrilateral vortex ring, control points on the end 
cap panels are close to the upper and lower surface panel edges for 
certain values of A0* This effect may lead to local numerical 
problems in the solution. It is clear that no pitching moment 
correlation can be obtained using this approach. 
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Figure 2.3: Semi-span wing and tip modeled as separate components 





Figure 2.4: Inboard and outboard sections showning 

tips 
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Figure 2.5: Tip lift and pitching moment variation with angle 

of attack; separate components, closed tips 
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The next step consisted of removing the panels on the wing and 
tip walls forming the gap, while maintaining the finite distance 
between the two surfaces. According to Reference 4, this type of 
approach is not allowed, since some "leakage" of the internal flow 
potential will occur. This will also prevent the boundary 
conditions from being satisfied over the entire surface and internal 
volume; but if the end caps are small, the effect on the integrated 
force and moment coefficients should be minimal. Figure 2.6 shows a 
substantial improvement in the prediction of lift with this new 
approach. A major improvement is also noted in the pitching moment 
curve, where the standard linear behavior is now apparent. However, 
the predicted theoretical slope is more positive than the 
experimental value. 

Since the integrated tip lift prediction is excellent using 


this approach, the data suggest an error in the prediction of 
the tip surface aerodynamic center location, X AC* Calculation of 


X Ar by taking the ratio of (C ) T to (Cr )t shows that the value 
obtained from VSAERO results falls short by 14.6% of the 


experimental value. VSAERO thus predicts X^q to be closer to the 

leading edge, thus resulting in a smaller value of (C m ) T than for 

a 

the experimental data. As will be seen later, this error in X AC 
prediction is not unique to method 2 (no end cap panels). It is 


difficult to determine the reason for errors in prediction. One 
possible cause may be related to the assumption that no separation 


takes place over the tip surface for the range of studied. If 
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Figure 2.6: Tip lift and pitching moment variation with angle 

of attack; separate components, open tips 
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separation did occur in the actual wind tunnel tests at higher *■ 
values of o&p the code could not predict with great accuracy the 
surface pressure in the separated region given the attached flow 
model used* Clearly, the improvements in and vs Oj. observed 
in Figure 2.6 as compared to the results of Figure 2.5 are related 
to the removal of panels from different components facing each other 
across the gap. Also, as will be shown in Section 2.2.3, the skin 
friction drag component is consistently underpredicted by the 
version of VSAERO available at NASA Ames. Some error thus results 
in the computation of C , since drag contributes to a small extent 
to pitching moment (see Appendix A, Equations A. 21 and A. 25). The 
limitations described here indicate that satisfactory prediction of 
both and using the same modeling technique is difficult to 
achieve. Despite this, the third modeling method was 
investigated. In this approach, the configuration consists of a 
single continuous surface. Figure 2.7 shows that the inboard wing 
and tip section are linked by a common section made up of panels 
which are stretched into a near vertical position. Figure 2.8 
illustrates the resulting and curves compared to the 
experimental data and the open-tips method results. 

The overprediction of lift can be attributed to the presence of 
flow singularities on the panels spanning the gap. Since a 
continuous wake was defined for the entire span in this modeling 
method, the Kutta condition was imposed at the trailing edges of the 
panels across the gap. Because of this condition, the solution 
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Figure 2.7: Semi-span wing and tip modeled as a continuous 

component with step change in incidence 
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Figure 2.8: Tip lift and pitching moment variation with angle 

of attack — TIP 1A 
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yielded doublet strengths in this area which were roughly of the 
same magnitude as those located on the adjacent lifting surfaces. 
These flow singularities induce velocities which alter the angle of 
attack distribution at the tip, thus resulting in an overprediction 
of lift. Although the resultant force normal to these panels 
produces some lift, its contribution was not added to the tip 
lift. The pitching moment curve is also steeper than that obtained 
from experimental results. This indicates that the overprediction 
of lift overcomes the error in X AC to yield higher magnitudes of 
pitching moment • 

Figure 2.9 shows the variation of tip pitching moment with tip 
lift. The underpredicted slope obtained from the open- tips method 
can again be attributed to the error in X AC prediction. As 
expected, this plot also confirms that VSAERO does not predict the 
zero-lift pitching moment accurately, thus resulting in a vertical 
offset of the curve compared to experimental results. Note the 
excellent slope correlation associated with the step change in 
incidence method. 

To verify the trends observed in lift and pitching moment 
variation for TIP 1A, a comparison of the characteristics of TIP IB 
(Figure 2.10) with experimental data was carried out. Notice that 
TIP IB has a higher sweep angle, and the aerodynamic center is 
farther aft than for TIP 1A. The variation of lift and pitching 
moment coefficients with angle of attack are shown in Figure 2.11. 
The same basic trends observed for TIP 1A also prevail in this 
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Figure 2.9: Tip pitching variation with lift — TIP 1A 
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Figure 2.10: Semi-span wing configuration with TIP IB 
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Figure 2.11: Tip lift and pitching moment variation with angle 

of attack — TIP IB 
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case. Although the lift prediction is more accurate with the open- 

tips method, it is slightly underpredicted. Because of the higher 

quarter chord sweep angle of this tip planform, it is likely that 

larger regions of separated flow will occur at high angles of 

attack, as compared to TIP 1A. This would result in higher levels 

of pressure drag. The underpredicted magnitudes of Cr along with 

U T 

the error in X AC location, poor skin friction and pressure drag 
predictions, probably all acount for the degradation of the pitching 
moment correlation observed here for the open-tips method. 

Once again, the step change in Incidence yields slightly higher 
lift and pitching moment variations with angle of attack. Figure 
2.12 also confirms the previous conclusion concerning the dC^/dC^ 
slope: the single-component method yields the best correlation. 

Note that since the. lift is overpredicted in the vs aj. curve 
(Figure 2.11), it is reasonable to assume that the good correlation 
shown in Figure 2.12 for the step change in incidence approach 
(l.e., almost matching slopes) is due to a combination of the 
overprediction of lift in the positive range (negative C m ^) and 
a resultant X AC closer to the pitch axis than would be obtained 
experimentally. As mentioned previously, the overprediction of lift 
is related, to the induced velocities at the tip. In spite of these 
problems, the data of Figures 2.8, 2.9, 2.11, and 2.12 show that the 
continuous component with step change in incidence yields the best 
prediction of pitching moment slope. For the range of sweep angles 
investigated, the correlations show that a maximum error of roughly 
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Figure 2.12: Tip pitching moment variation with lift — TIP IB 
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10-15% can be expected in dC^ 
value. 


/dctj. as compared to the experimental 


2.2.2 Spamrlae Loading Correlation 

To analyze the correlations of the previous sections in further 
detail, a study of the spanwise variation of sectional lift 
coefficient was carried out. Experimental data obtained from a wind 
tunnel test of the same configuration under consideration here were 
obtained from Reference 6. Flowfield velocities were measured using 
a Laser Doppler Velocimeter apparatus. The circulation at specified 
span stations was calculated by integration of these velocities 
along closed paths. The comparisons presented here are for TIP 1A. 

The VSAERO data were obtained at 12 spanwise stations for the 
Inboard wing and at 10 stations in the tip region. The numerically 
predicted values of sectional lift were normalized using the inboard 
wing chord. Figure 2.13 illustrates the correlation of these data 
with experimental values for A0 * 0.7°, and -4.3°. Results obtained 
from the three modeling methods discussed in the previous section 
are shown. The Inboard wing angle of attack is 6.6° for all four 
cases; the structural discontinuity is located at Y/(B/2) = 0.676. 
Figure 2.13 clearly shows the singular behavior at the structural 
discontinuity when the closed-tips method is used. As explained in 
Section 2.2.1, this is caused by potential flow effects and/or 
control point/panel edge proximity. For positive values of tip lift 
coefficient, it is seen that the step change in incidence method 


37 





38 






overpredicts the sectional lift, resulting in higher values of 
as shown in Figures 2*8 and 2.11. The sudden spanwise drop in 
for method 2 at Y/(B/2) = 0.68 is related to the structural 
discontinuity there. Since there is no surface to support a load 
and no singularity distribution over the gap, the lift must go to 
zero. These plots confirm that despite the marginal degree of 
correlation in the region near the gap, the open-tips method yields 
the best integrated tip lift coefficient for this lower range of A9. 

Figure 2.14 (A0 = -8.7° and -12.4°) shows that as the incidence 
angle of the tip relative to the inboard wing increases, the 
singularity peaks at the gap decrease in magnitude. In fact, it 
appears that for large values of A0, the closed-tips method would 
yield the best prediction of tip lift (Figure 2.14, A0 =* -12.4°). 
Overall, the step change in incidence method, selected for pitching 
moment prediction (Section 2.2), usually yields a more positive tip 
lift prediction than measured experimentally. 


2.2.3 Drag Correlation 

The purpose of this section is to evaluate the drag estimation 
capability of VSAERO. This is necessary to acquire an appreciation 
for drag related data and associated discussions presented 
throughout this report. 

In the present application of the program (assuming attached 
flow for all lift coefficients), the code predicts the drag due to 
lift and the skin friction drag. Unfortunately, the version of 
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Figure 2.14: Spanwise loading correlation, A0 = 8.7°, -12.4° 
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VSAERO available at NASA Aines during the spring of 1986 only printed 
out the skin friction drag as computed from the first iteration; 
this parameter was not updated in subsequent iterations. 

Figure 2.15 shows the free-tip drag polars from VSAERO data 
only, experimental data, and a combination of VSAERO drag due to 
lift and skin friction drag predicted empirically. The 
configuration was modeled here as a continuous component with a step 
change in incidence angle at the tip. Observe from Figure 2.15 that 
the drag calculated by VSAERO only is highly underpredicted. As 
discussed in Section 2.2.1, this is believed to have a small effect 
on the correlation of pitching moment. The empirical relation used 
is based on the flat plate analogy and is given by (Reference 7) as 

C„ - M 1 + L < t/C > + 100 < t/c >XsW S Rl!F (2 ' 2) 

SF 

For the particular airfoil section used, L = 2.0 (maximum thickness 
point is not located beyond x/c = 0.30), and R^g = 1*0 in 
incompressible flow. Fully developed turbulent flow was assumed. 
Figure 2.16 shows the same drag polars but with the configuration 
modeled as two separate components with open tips at the gap between 
the surfaces. Although this method seems to predict the lift 
characteristics with better accuracy, a comparison of Figures 2.15 
and 2.16 indicates only a slight improvement in the drag variation 
with lift correlation. This improvement is more apparent for the 
TIP IB polar at positive lift coefficients. 
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Re# = 0.867 million 


Figure 2.15: 


Drag polar comparisons — single component, step change 
in incidence, A0 =* -5° 
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Re# =■ 0.867 million 


Figure 2.16: 


Drag polar comparisons — separate components, 
open tips, A0 ■» -5° 
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Obviously, the version of VSAERO which was used does not have 
the capability to predict total drag accurately. It is seen that 
when the empirically obtained drag component is added to the induced 
drag, the minimum drag values for both tips are in good agreement. 
Also notice that the shape of the two curves is essentially the 
same, but the tip C^'s at which minimum drag occurs do not correlate 
well: VSAERO seems to overpredict the minimum drag in both 
cases. This could be related to interference effects which are not 
properly accounted for in the VSAERO model. Also, the empirical 
expression (Equation 2.2) is valid only for small incidence angles, 
thereby introducing more error at high angles of attack. 

For obvious reasons, it was decided at this point that only the 
drag resulting from the streamwise component of the pressure 
coefficient acting on the panel area (i.e. drag due to lift) should 
be used from VSAERO output. In subsequent computations of total 
drag, the skin friction component computed by the empirical 
relationship will be added to the non-skin friction contribution of 
VSAERO. 

As an illustration of this approach, consider the variation of 
total tip drag with lift coefficient squared. Figure 2.17 shows 
this relationship for both TIP 1A and TIP IB. In both cases, the 
tip at negative incidence actually produces less drag than the 
undeflected tip (A0 = 0.0). This confirms the hypothesis of 
Reference 2, which stated that the upwash caused by the vortex shed 
at the wing/tip junction reduces the streamwise component of normal 
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Figure 2.17: Variation of tip drag with lift — TIP 1A and TIP IB 
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force at spanwise stations closest to the tip vortex, thereby 
reducing the total induced drag of the tip. This is one of the 
beneficial features of the free-tip concept and should result in a 
drop in the rotary wing power requirement. Note also that the 
predicted zero-lift for both AOs is the same, further supporting 
this theory. 

2.3 VELOCITY VECTOR PLOTS 

As shown by Equation (A. 18) in Appendix A, the velocity at any 
point in the flow field may be computed by the VSAERO program. The 
user specifies points in 3-D space at which the velocity and 
pressure coefficient are to be calculated. With the three velocity 
components known at the specified points, a plotting routine can be 
used to obtain vector plots showing the projection of the resultant 
total velocity in a given plane. This approach can be used to 
visualize the development of tip vortices and determine the relative 
magnitude of the velocities induced by vortex flows. As will be 
shown in Chapter 3, the velocities induced by three-dimensional 
effects contribute to drag reduction in certain cases. Vector plot 
results could also be used in conjunction with experimental 
techniques for visualization of flow conditions. 

The method consists of specifying "Velocity Scan Planes" which 
intersect the two surfaces under consideration and are perpendicular 
to the free stream velocity direction. In this case, eight planes 
were selected to analyze the flow over TIP 1A. The first plane is 
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located immediately ahead of the leading edge while the last is 
roughly one-third of a chord length behind the trailing edge. 

Figure 2.18 shows the location of the eight equally spaced planes. 
Note that they overlap the structural discontinuity of the semi-span 
wing and the outboard tip, which are the primary regions of 
interest. A total of 400 grid points were generated. Figure 2.18 
also depicts the two cases which will be discussed. 

Figure 2.19 shows the projected velocity vectors in the y-z 
plane for = 12.0° and = 2.0°. The structural discontinuity is 
located at y = 65.02 cm, and the outboard tip at y = 95.77 cm. From 
the initial perturbations at the leading edge to the region aft of 
the trailing edge, the velocity vectors calculated from the 
distribution of singularities on the body and wake illustrate the 
development of the inboard and outboard vortices. Although the 
lifting surfaces are interacting, two distinct tip vortices are shed 
at the structural discontinuity; this has been confirmed by flow 
visualization studies performed during the January 1985 test. With 
both surfaces generating positive lift, the circulation in the tip 
vortices of both surfaces will be of opposite sign. The interaction 
of these two vortices depends on the respective loading of each 
surface. For example, when 3 12.0° and aj» = 2.0°, the inboard 
wing tip vortex dominates the flow field around and aft of the gap 
(Figures 2.19, planes if 6-//8). However, its effect on the outboard 
surface is reduced slightly due to the vortex of opposite strength 
shed by that surface. 
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In the second case (Figure 2.20), the outboard tip section is 
producing negative lift, thus accounting for the clockwise direction 
of the circulation in the most outboard vortex (at y » 95.8 cm). At 
the structural discontinuity though, the vortices will merge 
together, since the circulation in both of these is of the same 
sign. This is confirmed by the large magnitudes of the induced 
velocities at the lower surface of the outboard section (Figure 
2.20, planes Z/4-//6, to the right of y a 65.02). These induced 
effects will be discussed more thoroughly in Section 3.3.1. 
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Figure 2.20: Vector plots for TIP 1A, cfy = 5°, 


5° (continued) 
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3. ANALYSIS OF NEW TIP PLANFORMS 


The major objective of this phase of the research is to design 
and analyze new tip planform shapes for the free-tip. The emphasis 
in this process is on moving the aerodynamic center aft to improve 
the response of the free-tip to angle-of-attack variations. These 
new tip shapes are all related through a common inboard half. The 
quarter chord sweep angles as well as the taper ratios of the tips 
were varied, leading to the desired aerodynamic effects. The free- 
tip performance characteristics to be investigated are: 

© pitching moment variation with angle of attack 

0 compressibility effects on and for representative 
subcritical Mach numbers 

© sectional and surface aerodynamic center shift at high 
subsonic Mach numbers 

# tip planform drag and its effect on the semi-span wing 
drag polars. 

In this chapter, the computed values of are normalized by 
the inboard wing chord ( C XNB ) * Chapter 2, the moment 

reference point is the tip pitch axis, which coincides with the 
inboard wing quarter-chord line. 

3.1 PITCHING MOMENT VARIATION OF THE NEW PLANFORMS 

Table 3.1 lists the basic geometric characteristics of the new 
tip planforms. The details of the geometry are illustrated in 
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Table 3.1: Geometric Parameters of the New Tip Planforms 
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Figure 3*1 • The new tips are compared here to the old baseline 

configuration described in Chapter 2 (TIP 1A). A new baseline 

configuration [TIP 2, Figure 3.1(a)], which establishes the 

foundation on which the new planforms are based, was designed 

strictly on the requirement to Increase the slope of the vs oj, 

curve* Since the inertia properties of the free-tip and the torque 

output of the passive controller are unknown at this time (these are 

some of the parameters along with (C m ) T which appear in the 

a 

equations of motion of the tip, Reference 3), tips with intermediate 

values of (C m ) T were designed. These might be useful to optimize 
a 

future designs once all the relevant response characteristics are 
known. The incompressible pitching moment behavior of the new tips 
is compared to that of the old baseline configuration (TIP 1A) in 
Figure 3.2. 

As may be observed, increasing the sweep angles and moving the 

aerodynamic center aft has clearly resulted in a negative increase 

in (C m )f. This obviously Improves the response of the new tips to 
a 

variations in local angle of attack or inflow. Table 3.2 lists the 
magnitudes of (C m ) T for the tips of Figure 3.2. 
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Figure 3.2: Comparison of pitching moment variation with angle of attack for new tip planforms 



Table 3.2: Variation of Pitching Moment with Angle of Attack 

for New Free-Tip Planforms 




( C m >T 
a 


Tip // 

(deg -1 ) 

1A 

(old baseline) 

-0.0073 

2 

(new baseline) 

-0.0165 

3 


-0.0145 

4 


-0.0125 

5 


-0.0115 

6 


-0.0111 


It is seen that from the old baseline configuration to TIP 2, a 

125% increase in (C m ) T has been achieved. Also, tips 3 through 6 

a 

provide a good range of intermediate values of this parameter. As 

stated earlier, these could be used for future designs. Because of 

the small difference in (C m ) T between tips 5 and 6, the former will 

a 

be dropped from the analyses that follow in this report. 

3.2 COMPRESSIBILITY EFFECTS 

3.2.1 Prediction of Tip Lift and Pitching Moment Data 

The installation of free-tips on the rotor blades of a 
helicopter should result in power and vibration reductions in 
forward flight. The highest tip Mach numbers are encountered at 
high advance ratios. It is therefore necessary to investigate the 
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performance of the decoupled tip of a semi-span wing at high Mach 
numbers, since prevailing compressibility effects will alter the 
aerodynamic characteristics of the tip. Unfortunately, the VSAERO 
program is a subsonic code and may therefore be used only for a 
combination of onset velocities and lift coefficients which do not 
produce shocks over the surface of the tip. 

If the free-stream Mach number were increased above this range, 
shocks could form on the body upon local deceleration to subsonic 
flow. When the flow becomes supersonic locally and normal shocks 
form, flow separation is primarily due to the large adverse (and 
steep) pressure gradients associated with these shocks, and Reynolds 
number effects become secondary. VSAERO cannot predict this form of 
separation phenomena, since it cannot predict the formation of the 
shock itself; therefore, it is desirable to avoid these critical 
onset conditions. Note that for free-stream Mach numbers above the 
shock-free limit, VSAERO still yields continuous chordwise pressure 
distributions; but because flow separation is not predicted, the 
lift coefficients resulting from chordwise integration of Cp are 
overpredicted. 

The method used in the code to account for compressibility 
effects is straightf orward. After the pressure coefficient at the 
panel center has been determined for incompressible flow (see 
Appendix A, Equation A. 22), it is corrected using a simple 
compressibility factor. For the cases discussed in this chapter, 
the Karman-Tsien option in the program was used. This relation is 
given by 
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Note that for low M^, this expression is equivalent to the Prandtl- 
Glauert factor. Once the Cp at each panel has been corrected for 
compressibility, the chordwise and spanwise summations are carried 
out in the manner described in Appendix A. 

In the study of shock-free compressibility effects, a realistic 
range of Mach number and tip needs to be investigated. Figure 
3.3, taken from Reference 2, illustrates the range of A0 as a 
function of azimuth angle obtained from the first experimental test 
of a rotor equipped with free-tips. The largest contrast in 
aerodynamic effects for a helicopter in forward flight is found by 
comparing the ij; = 90° (advancing side, highest tip Mach number) with 
the « 270° (retreating side, lowest tip Mach number). For this 
analysis, it will be assumed that at = 90°, the pitch angle of the 
inboard blade at the structural discontinuity (R = 0.9) varies 
between 0° and 4°. At the = 270° position, the local pitch angle 
is assumed to vary between 10° and 12°. These are representative of 
the blade pitch at this radial location for typical advance ratios 
(Reference 8). Therefore, if these incidence angles are used for 
the semi-span- wing inboard section, the equivalent tip angles of 
attack may be obtained from the A0 values of Figure 3.3. 

As discussed earlier, if a subsonic code such as VSAERO is used 
in this analysis, it is necessary to remain below free-stream Mach 
numbers and tip C^ f s which lead to the formation of shocks. Two 





sources of data were used to ensure that compression shocks did not 
occur. The first consisted of checking the local Mach number on the 
upper and lower surfaces of the tip as predicted by VSAERO. The 
code computes the Mach number at panel centers as a function of 
corrected Cp values (for compressibility) and free-stream Mach 
number. The second involved the use of 2-D experimental data 
(Reference 9). For specified values of M^ and C^, these data were 
used to avoid lift coefficients in the nonlinear range. In an 
attempt to compare the differences between typical advancing and 
retreating side values of tip Mach number for a small-scale rotor, 
two cases were considered, M^ = 0.49 and M^ = 0.70. The method used 
to generate the values of and is described in Appendix B. 

After several attempts in identifying the highest possible 

combinations of M_ and C T , it was found that the critical case 
°° L.^ 

occurred at M^ = 0.49 on the inboard portion of the semi-span 
wing. At = 0.49, there still is a good safety margin away from 
nonlinear effects for the tip region. The final combination of 
parameters which yield the shock-free flow over the entire semi-span 
wing at M^ = 0.70 and M a = 0.49 are listed in Table 3.3. 


Table 3.3: Combinations of Pitch Angles for Shock-Free Flow 

at the Selected Free-Stream Onset Velocities 


Advancing side: M^ = 0.70 

0.0 < cfo < 4.0, oj = -8.0°, -5.0°, -2.0° 
Retreating side: M^ = 0.49 

8.0 < < 12.0, ouj, = -14.0°, -11.0°, -8.0° 
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Note that the test RPM and geometry of the rotor used in the 
tests of Reference 2 were used to evaluate the related Mach numbers 
in the tip region on the advancing and retreating sides. Tip Mach 
numbers of 0.70 and 0.49 correspond to an advance ratio of y - 0.176 
for the small-scale rotor in the wind tunnel. Since only the 
critical cases (with respect to VSAERO capabilities) were of 
interest in this study, Mach numbers which would correspond to lower 
values of y were not considered. Also, it was found later that the 
limits of Table 3.3 did not change with respect to the various tip 
planforms considered; therefore, these incidence angles are valid 
for all of the planforms of Figure 3.1 at the pre-selected limit 
Mach numbers (M^ = 0.70, 0.49). 

Despite the precautions listed in the previous paragraphs, flow 
separation can still occur due to Reynolds number effects. Since 
previous attempts to stitch the wake on the lower surface of the tip 
failed, the separation line was specified at the trailing edge for 
all cases. If values of and in the nonlinear range (i.e. 
with partially separated flow) had been sought, an exhaustive 
iterative procedure would have been initiated. This would have 
required accurate specification of the predicted separation line, 
based on output from the boundary layer routine. The data completed 
for cases where separated flow is observed thus result from a first 
estimate only. These were not obtained from converged solutions, 
and are presented to verify trends only. The predicted values of 
C T and (L, for the conditions listed in Table 3.3 are tabulated in 

Uip Ulm 

Appendix C. 
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Table 3.4 lists the low angle of attack lift and pitching 

moment slopes for tips 2, 3, 4, and 6 calculated from the data in 

Appendix C. Only values of and which fell in the linear 

range were used to estimate C L and C m of the tip. These are 

a a 

compared to the computed values in incompressible flow. The 

following observations can be made: 

First, note that the incompressible lift curve slopes (Or ) T 

a 


are higher than would be expected for such a low aspect ratio 

lifting surface (A T » 1.80). This is consistent with the results 

obtained for the step change in incidence method (see Section 2.2.1) 

which was used here to obtain better (C m )^ predictions. Increases 

a 

in (C L ) T more closely agree with the Prandtl-Glauert factor for 
a 

= 0.49. The larger deviation at M ot = 0.70 is expected, since 

experimental data show that the Prandtl-Glauert law may overpredict 

C L at higher subsonic Mah numbers (Reference 10). 
a 

Comparison of the % increase in (C )rj, from M w s 0 to 0.49 and 

a 

from = 0.49 to 0.7 shows that compressibility effects are more 

pronounced on the pitching moment slope (C m ) T at = 0.70. This 

a 

can be partially attributed to the aft shift in aerodynamic center 

with Mach number. When the total moment is computed, the resultant 

lift acting at the A.C. has increased due to compressibility, but so 

has the moment arm to the reference point (the shaft pitch axis in 

this case). However, at = 0.49, (C m does not increase as 

a 

much. In fact, the increase closely follows that of the tip lift 
curve slope, (C L )^, indicating that no shift in the A.C. has 
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Table 3.4: Compressible Lift and Pitching Moment Slopes (Linear Range) 
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occurred. This is confirmed by the A.C. locations on the right hand 
column of Table 3.4 calculated from the lift and pitching moment 
slopes. When M a * 0.50, the position of the A.C. is still very 
close to its incompressible location. As will be shown in the next 
section, the sectional A.C. moves aft somewhere between = 0.5 and 
= 0.7 . 

00 

Because of the difference in (C m a large difference in tip 

CL 

response between the \|> = 90° and = 270° azimuthal positions can be 
expected. The average Increase for the four tips is 26% at 
p = 0.176, for a small-scale rotor (Reference 2). 

No particular tip seems to have a larger increase in (C m ) T 

Ot 

(compared to the incompressible case) than another at the same Mach 

number. The percentages in increase of (C m ) T are all within about 

a 

5% of each other. Any variations in the increase of this parameter 
may be related to uncertainties associated with slope estimation. 

3.2.2 Sectional Aerodynamic Center Shift 

The sectional aerodynamic center location is computed by using 

a relation derived from the sum of the moments about an arbitrary 

point on the section chord. With the reference point at the quarter 

chord of the section and since by definition, dC /da = 0, the 

m ac 

following expression can be obtained: 

x ac /c - 0.25 - dC m /dC t | 0 . 25c (3.2) 
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Since the VSAERO program computes the sectional lift and 
pitching moment about the quarter chord at each user-defined 
spanwise column, the slope dC m /dC ^ | q . 25c can be easil y obtained for 
a specified Mach number. The results for Mach numbers up to M a = 
0.80 are shown in Figure 3.4. Using TIP 6 as an example, the 

computations were carried out at the two span stations indicated at 

the top of the figure. To maintain shock-free flow, only values of 
sectional lift coefficient between 0.10 and 0.25 were used to 
compute x ac /c. Although the flows here are highly three- 
dimensional, qualitative comparisons with 2-D experimental results 
can be made. The data of Figure 3.5, taken from Reference 9, 
illustrate the shift in aerodynamic center as a function of free 
stream Mach number for the V23010 airfoil. At both span locations, 
the trends of Figure 3.4 are similar to those observed in Figure 
3.5. The forward travel of the A.C. around M^ = 0.60 to 0.65 is not 

as significant as observed from the experimental data (« 3.4% of the 

chord); however, the Mach number at which the A.C. begins its aft 
shift matches the 2-D value quite well. Due to the limitations of 
the theoretical model, which is incapable of predicting shocks and 
their effects on pressure distribution, the data predicted by VSAERO 
are unreliable beyond M w = 0.72. 

Figure 3.4 also confirms that the computed results are in 
agreement with the theory of finite aspect ratio wings. The 
following relation, taken from Reference II, yields the sectional 
aerodynamic center location for an assumed wing loading function: 
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Figure 3- 



SECTIONAL AERODYNAMIC CENTER SHIFT 
WTTH FREE STREAM MACH NUMBER 



: Compressibility effects on sectional aerodynamic center 

location, V23010 airfoil — VSAERO data 
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Figure 3.5: Compressibility effects on aerodynamic center location 

V23010 airfoil — experimental data taken from 
Reference 9 





(3.3) 


x ac /c = 1/2 “ n > 

The parameter "n" is a function of aspect ratio: the value n = 1/2 

represents the limit A > and the limit n -► 1 represents the case 
where A « 0. As the data from Figure 3.4 for = 0 indicate, the 
sectional A.C. of the inboard wing is in the vicinity of the quarter 
chord. Because symmetry about the wing root was assumed in this 
case, the inboard portion behaves as a high aspect ratio wing, and 
the predicted A.C. location in incompressible flow is in good 
agreement with Equation (3.3). The theory also predicts that for 
low aspect ratio wings, the sectional A.C. moves towards the leading 
edge. The computed data at = 0 confirm that this occurs here 
also, because of the three-dimensionality of the flow in the tip 
region. 

Figure 3.6 illustrates sectional chordwise pressure 
distributions at M a = 0.0, 0.6, and 0.7. Note that these Cp 
distributions are plotted for the same sectional lift coefficient; 
the span stations correspond to those where the aerodynamic center 
locations were evaluated. The shifts in aerodynamic center location 
observed in Figure 3.4 between = 0.0 and = 0.70 are of the 
order of 1-2% of the chord length. This finding is correlated by 
the data of Figure 3.5 at = 0.60 and 0.70, where the peak Cp 
values have moved aft slightly. This can be observed also by 
comparing the steep negative increases in Cp at the leading edge. 
These portions of the curves for = 0.60 and 0.70 are aft of the 
solid line (incompressible flow). Since the sectional C^ is 
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I 

Figure 3.6: Compressibility effects on sectional pressure 

distribution, V23010 airfoil— VSAERO data 
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constant for all three Mach numbers, this aft shift in the loading 

(along with an increase in C ) accounts for the negative increase 

ac 

in pitching moment. Also, notice that as is increased, the 
gradients in C p are steeper aft of the upper surface pressure 
peaks. This is where shocks would most likely occur if the 
sectional loading were increased. The small peaks in C p on the 
lower surface at the leading edge are probably associated with the 
drooped nose of the V23010 airfoil. The undeflected trailing edge 
tab produces the flat C p distribution over the last 10% of the 
chord. 


The position of the A.C. of a three-dimensional surface is 
defined as the weighted average of the sectional aerodynamic center 
locations as follows: 


X 


AC 


b/2 

/ 

-b/2 

~rrr 


c A cx dy 
l ac J 
a 


(3.4) 


-b/2 *a 

where C 9 is the sectional lift coefficient due to the local angle 
*a 

of attack, c is the local chord length, and b is the lifting surface 
span. In Table 3.4, the surface A.C. locations as a fraction of the 
inboard wing chord are listed. By comparing the 3-D shifts between 
= 0.5 and 0.7, it can be seen that the A.C. travels aft by 1.7 to 
3.2% of the inboard wing chord length for the four tips shown. This 
seems to be reasonable in view of the calculated sectional AC shift 
for TIP 6. 
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The data presented in this section were obtained through the 

use of the maximum capabilities of the VSAERO code in compressible 

flow analysis. Because of the limitations of subsonic lifting 

surface theory and the approximation of compressibility effects, 

some uncertainty is associated with these results. In addition, 

many of the predicted parameters only have limited ranges of 

application. For example, the (C^ and (C m ) T slopes of Table 3.4 

ot a 

are valid only for a small range of tip angles of attack or lift 
coefficients. Nonetheless, the computed data show good consistency 
with sectional and three-dimensional compressible flow theory. 

These results can be taken into account in the design process of tip 
planforms for the free-tip rotor. 


3.3 DRAG POLARS OF NEW TIP DESIGNS 

The drag of a semi-span wing with a structurally decoupled tip 
is important in this analysis because of its relationship to the 
power requirement of a helicopter rotor. In rotary wing vehicles, 
the power requirement is directly related to the torque created by 
the drag of the blades. The tip region of the blade is of primary 
interest because of its longer moment arm to the helicopter hub and 
associated high dynamic pressures which Increase the drag 
considerably. In this section, the drag polars of two 
representative tip planforms will be presented and analyzed. The 
effect of tip drag on the complete semi-span configuration will also 
be investigated. 
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The configuration is modeled as a single component with a step 
change in incidence at the tip. The skin friction drag is estimated 
empirically and added to the induced drag predicted by VSAERO. The 
reader is once again reminded of the main assumption in this 
analysis, in that despite the constant spanwise dynamic pressure 
limitation of VSAERO, there exists a relationship between the data 
produced for a semi-span wing and the performance of a free-tip on a 
helicopter rotor blade in forward flight. 

3.3.1 Tip Planform Drag Polars 

The drag polars for TIP 2 and TIP 6 are shown in Figure 3.7. 

The tip planform areas were used to normalize the calculated lift 
and drag. Fewer data points were computed for A0 ■ +5.0° to remain 
within the linear range of the lift curve. 

The data are presented for three values of A 9 to illustrate the 
beneficial aspects of the free-tip concept in incompressible flow. 
Clearly, the drag produced in the tip region is substantially 
reduced for negative tip incidence angles. During the rotor tests 
of Reference 2, power reductions were recorded for tip Mach numbers 
up to Hj « 0.85 for A0 < 0°. Due to the limitations in compressible 
drag prediction, this cannot be verified for the semi-span wing and 
tip configuration with the VSAERO theoretical model. 
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However, in incompressible flow, the drag reduction effect may 
be explained by considering a 2-D section cut of the tip surface. 

As shown in Figure 3.8, when the tip is undeflected, the downwash 
induced at the quarter chord by the semi-span wing tip vortex 
reduces the effective angle of attack of the section under 
consideration. The free-stream component of the resultant force 
(parallel to the onset velocity V^) is the induced drag contribution 
of the section. Now consider the case where the tip has an 
incidence of -5 degrees relative to the inboard section. A strong 
vortex is shed at the structural discontinuity due to the large 
spanwise lift gradient there. For a positive inboard loading which 
yields exactly the same sectional lift as for the A0 *» 0.0 case, the 
velocities induced by the inboard vortex will reduce the magnitude 
of the downwash at the section. As illustrated in the bottom 
portion of Figure 3.8, this will result in a lower Aa and hence, 
less induced drag produced at the section under consideration. When 
the drag components of each section are integrated, the induced drag 
of the tip at A0 = -5.0° is below the computed value for A0 = 0.0°. 

The effects illustrated in Figure 3.8 thus verify the tip drag 
reductions measured during the first semi-span test (Reference 1) as 
reported in Reference 2. The semi-span wing and flow conditions 
modeled in the present study differ from those of the experimental 
rotor configuration used previously (Reference 2); also, as 
mentioned earlier, compressibility effects are not taken into 
account in the present drag prediction process. Despite this, 
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the induced drag reduction effect is believed to contribute 

significantly to the power reductions observed during the tests of 

Reference 2, even at high tip Mach numbers. 

Since the analyses of TIP 2 and TIP 6 yielded equivalent drag 

polars (almost point for point), it was not judged necessary to 

produce similar data for TIP 3 and TIP 4 which are geometrically 

similar. This also shows that the tip geometry may be tailored to 

given performance requirements (i.e. a target value of C ) without 

a 

affecting the trends observed in the drag polar of the new baseline 
configuration (TIP 2). 


3.3.2 Semi -Span Wing Drag Polars 


Figure 3.9 shows the semi-span drag polars with different tip 
planforms as predicted by the method outlined previously (VSAERO and 
empirical data). The first effect which is readily apparent is that 
the drag polars are all identical; this follows from the results of 
Section 3.3.1, which showed little difference in drag behavior 
between the various tip planforms. Also, the tips contribute a 
fraction of the drag, due to the normalization of the integrated 
drag force by the total reference area of the semi-span wing. A 
slightly higher ratio of S^ET^REF ( i#e# higher skin friction drag) 
accounts for the drag increase of TIP 6 (Figure 3.9). 

The wing drag polars are plotted separately in Figure 3.10 for 
each tip planform with A6 = -5°, 0°, +5 9 . These illustrate a point 


of major interest when one considers the differences in 



with tip 
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R eif = 0.867 million 


Figure 3. 


: Semi-span wing drag polars with new tips at 

A0 = 0°, -5° 
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Figure 3.10: Semi-span wing drag polars with new tips 

(continues) 
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Re// - 0.867 million 


Figure 3.10: Semi-span wing drag polars with new tips (concluded) 
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incidence angle. Observe that at low wing C^'s, the undeflected tip 
configuration produces less drag than with the tip at negative 
incidence angles. Since it was shown in Section 3.3.1 that negative 
tip deflections yield less tip drag than zero-incidence 
configurations, it is apparent that a sharp increase in the semi- 
span's drag occurs for negative values of A0. 

These observations, although difficult to quantify, may be 

described qualitatively. When the tip is undeflected, the effective 

aspect ratio of the wing is greater than at A0 = -5.0; therefore, 

the induced drag of the semi-span wing at A0 = 0° will be at a 

minimum. When the tip has a finite deflection angle, the vortex is 

shed at the tip junction (increasing the induced drag of the wing) 

rather than at the outboard extremity. From Figure 3.10, for 

negative A0's at low values of (and thus, low Cj^) the reduction 

in tip drag is not sufficient to overcome the increase in of the 

wing. However, note from Figure 3.7 that at higher C, , the gap 

T 

between tip C^'s Increases (at constant C^) . Therefore, a point is 
reached where the drag reduction in the tip region balances the 
increase in wing drag; this is illustrated in Figure 3.10, where the 
curves for A0 = -5.0° and 0.0° merge. 

Since this approach does not account for the radial variation 
of dynamic pressure, one can only speculate about the resulting drag 
of a free-tip equipped helicopter rotor. In numerical analyses 
which use blade element theory, the sectional drag coefficients vary 
radially, as they are a function of local angle of attack, Reynolds 
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number, and Mach number. Whatever the local value of C^, the drag 
force computed for each blade element will increase in the outboard 
direction because of increasing dynamic pressure. Since the free- 
tip is located in the region of highest dynamic pressure, the drag 
reduction which takes place there for negative incidence angles is 
probably significant. This effect is the most probable cause for 
the drop in required power which was observed during the small-scale 
rotor test discussed in Reference 2. 
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4. CONCLUSIONS 


The following conclusions can be made regarding the analysis of 
the semi-span wing with VSAERO: 

• The VSAERO computer program has been used to model the 
characteristics of a semi-span wing with a deflected 
tip. The knowledge gained in using the code indicates 
that the best modeling techniques are as follows: 


The open-tips at the gap method gives best C^ vs 


°T 


correlation. 


• The continuous component with step change in 


incidence yields best C^ vs cij. and C^ vs C^ 
correlation. 


• The correlation of spanwise loading with experimental 
data is very satisfactory. 

o Drag prediction is unacceptable, primarily due to poor 
skin friction drag estimation. 

• The variation of C^ with odj may be increased considerably 
with modifications to the tip planform geometry, as shown 
by the design of new tips. 

• In shock-free compressible flow at = 0.70, the VSAERO 

model predicts an average of 47% increase in (C m from 

a 

its incompressible value. None of the new tips stand out 

with larger increases in (C m ) T at a given Mach number. 

o 

The predicted aerodynamic center travel is consistent with 
sectional and 3-D compressible flow theory. 
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a VSAERO analyses of a semi-span wing predict an increase in 

(C m of 26% when is increased from 0.49 to 0.70. 
a 

These Mach numbers are representative of retreating and 
advancing side tip Mach numbers for the small-scale rotor 
of Reference 2. 

• For negative tip incidence angles, it has been shown 
analytically that the tip drag is reduced. This confirms 
the experimental results of Reference 1. 

• At higher wing C^'s, the tip drag reduction balances the 
increase in induced drag of the inboard blade. 

• The drag polars of the semi-span configuration are not 
sensitive to planform geometry. 

Recommendation : 

For the next phase in the numerical analysis of the free-tip 
concept, consideration should be given to reproducing the exact 
conditions under which the free-tip should operate. This could be 
done by incorporating the equations of motion of the tip (including 
inertia effects, controller moment, and aerodynamic moment) to an 
existing rotor code which would account for the following effects: 

• Transonic Mach numbers (up to M^ = 0.9) 

• Unsteady aerodynamics and correction for viscous effects 

• Free wake capability with blade vortex interaction 
modeling. 

If a rotary wing code were used, a more accurate prediction of 
inboard blade effects induced on the tip would be possible. 
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APPENDIX A: 


VSAERO THEORETICAL FORMULATION 


The theoretical formulation of the code is based on the 
computation of the velocity potential at a point P in the flow field 
under consideration. The governing equations are (Reference 4) 

7 % = 0 

= 0 


*P 


1 _ 

4tt 


i i 

s+w+s 

00 


($ - $ i )n 


7 (i)ds 


- j ] ; n • (Vt - TO )dS (A. 1) 

S+W+S 

oo 

Where $ and are velocity potentials which exist outside and 
inside the region of interest, respectively, and "r" is the distance 
from the point P to the element dS on the surface; "n" is the unit 
normal vector to the surface pointing into the fluid (Figure A.l). 
The first integral represents the disturbance potential from a 
surface distribution of doublets with density ($ - per unit area 
while the second integral represents the contribution from a surface 
distribution of sources with density -n • (7$ - TO^) per unit 
area. The doublet density accounts for the local jump in potential, 
and the source density is associated with the local jump in normal 
components of velocity. 


A-l 




Since onset flow conditions are assumed to exist at S , the 

oo y 

integrals taken over this surface will reduce to <j> , the velocity 

00 

p 

potential at point P due to the onset flow. Also, the upper and 

lower wake surfaces are assumed to coincide into a thin wake, so the 

term for that integral disappears. If entrainement into the wake 

surface is neglected, the jump in normal component n • (V$ - ) 

U L 

vanishes so the source terms associated with the wake also 
disappear. Equation (A.l) above then becomes 

* P - ^ 1 / (* - V n * 7( 7 )ds 

s 

- ^ i i 1 n * (V$ - V* i )dS 

s 

•+ ^ J J (*u - V n * v( i )dW + (A,2) 

If point P lies on the surface, the integral becomes singular 
there and the point must be avoided. This is easily achieved by 
locally modifying the path of integration around point P in the form 
of a hemisphere. In the limit as the sphere radius goes to zero, 
the local contribution to velocity potential is l/2($ - $i)p> which 
is half the local jump in potential across the surface at P. The 
first potential in the expression being on the side of the surface 
on which P is located, if P is on the inside surface then Equation 
(A.2) becomes 


A-3 



$ p = j i <* - v n • v 4 )ds - 1 <* - $ i)p 

S— P 

- ~ I I - # • (V# - V$,)dS 
4 tt j j r l 

d 

+ 4w i i ($ U " V n * v( 7 )dw + *» p (A ’ 3) 

The external Neumann boundary condition is applied through the 
following relation: 

n • V<fr » -V_ T on S (A. 4) 

N 

where V N is the resultant normal component of velocity relative to 
the surface. Nonzero values of V N are used to model boundary layer 
displacement effects . 

The wake surface W does not support a load; the local vorticity 

vector (for an element dW) associated with the doublet 

distribution ($ - $ ) is given by 

U L 

Y = -n x 7($ u - $ L ) (A. 5) 

From the Kutta-Joukowski law, the elementary force on dW, in the 
presence of the mean velocity V, is 

6F = pV x ydW (A. 6) 

Since $ x y = 0 must hold for the force to be zero, then from 
Equation (A. 5) 

V x {n x V($ - $ )} =0 (A. 7) 

U li 

after expansion, 

nV • V($ - 4> ) - V • nV($ - $ ) = 0 (A. 8) 

U L U L 


A-4 



This equation is satisfied when V • n = 0 , i.e. if the surface W is 

aligned with the local flow direction. Also, V • V($> - <fr ) =0 

must be satisfied, so that the gradient of the doublet distribution 

is zero along mean streamlines in the wake surface; therefore, the 

wake doublet distribution is constant along mean streamlines in the 

wake and is determined by the condition at the point where the 

streamline leaves the surface S. Throughout the iterative process, 

the Kutta condition is implied by shedding the trailing edge 

potential jump (<I> - $ ) as a constant down each streamwise line on 

U Li 

an initially prescribed wake surface. Once the solution has 

converged, the upstream edge of the wake (hence the trailing edge) 

will carry no load, satisfying the Kutta condition. 

The application of an internal boundary condition is what makes 

the VSAERO approach accurate, convenient to use, and robust 

(relatively insensitive to bad paneling schemes). The essence of 

the procedure is to ensure that when passing through the boundary S, 

the jump in potential from the internal flow to the external flow 

should be small. This requires a minimum of perturbation from the 

singularities and is achieved by specifying that the internal flow 

potential be equal to the onset flow potential. Writing the total 

potential as the sum of the onset and perturbation 

potentials $ = <f> + <f>, and applying the internal Dirichlet boundary 

00 

condition $.=<(>, Equation (A. 3) becomes (with $ - d> ) 

« i OO P T oo 
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(A.9) 


J J < *« • V(i)dS - j <f> p + j / (* D - * L )m • V(i)dW 

S "”P w 

- J j - u • (V4 - d s = 0 

s r 

Comparing the above expression with Equation (A. 3), it is 
observed that now the doublet density may be expressed as 

4ttu * <J> = $ - <p m (A. 10) 

while the source distribution is 

4tto = - n * (V$ - V<f> ) (A. 11) 

00 

The first term in Equation (A. 11) may be replaced by the Nejrinann 

boundary condition (Equation 2.4) to obtain 

4tto - V - n • V (A. 12) 

N 00 

For the first potential flow solution, no boundary layer 
parameters are available; and since the boundary is solid, = 0. 
With the panel geometry already established, the terms n • V may be 
computed easily so the source distribution is known at the outset. 
Once the first potential solution is determined, the pressure 
distributions are transferred to the boundary layer routine. In 
subsequent potential flow calculations, the source distribution is 
given by 

- h ( V*> + V K0RM - » • ». < 4 ' 13 > 

where the first term on the RHS is obtained from boundary layer 
calculations • 
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The integral relation (Equation A. 9) may be solved for the 
unknown velocity potential distribution over the body surface. In 
the numerical procedure, the equation is satisfied at a finite 
number of points on the surface (one control point on each body 
panel with the doublet and source distributions assumed constant on 
each panel). The surface integrals therefore become summations over 
all panels, and the integral relation (Equation A. 9) is transformed 
into a set of simultaneous linear algebraic equations where the 
unknowns are the doublet strength on each body panel (recall that 
the source strengths are known — Equations 2.12 and 2.13): 


N N N 

j/w + + ' 01 J ' l * NS 

K=1 L=1 -L K — 1 


(A. 14) 


At the outset, the wake panel doublet strengths Uy are not known 

Ii 

and the second summation is not required. For subsequent 
iterations, the computed trailing edge potential jump py = 

Is shed as a constant down each wake column, as discussed before. 

In the code, the system of equations above is solved by a direct 
method for a moderate number of unknowns and by a blocked Gauss- 
Seidel iterative procedure for larger numbers of panels. Note that 
since <j> = Attp^, C = -2tt accounts for the — 1 / 2 <f> p term when the 
point P is on the surface. Ng and Ny are the number of surface and 
wake panels, respectively. Bj^ and Cj^ are the perturbation 
velocity potential influence coefficients for the constant source 
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and doublet distribution of panel K acting on the control point of 
panel J. From Equations (A. 9) and (A. 14) it is seen that 


B 


JK 


and 


1 J i dS 

Panel K 


(A. 15) 


= J j n • V(i)dS (A. 16) 

JK Panel K r 

The vectors n = n^ are established by the local geometry and are 
constant for each panel; r is the length of the vector from the 
surface elements of dS of panel K to the control point of panel J. 
When the doublet solution is known, Equation (A. 3) can be rewritten 
as 


$ = J J V n • V(i)dS + Ky + j j f dS 

S S 

+ J JUrj n • V(i)dW + (A. 17) 

W P 

where \i and a are defined by Equations (A. 10) and (A. 12) 
and \iy = are wake doublet values (constant for each wake 

column). If point P is off the surface, K = 0; if P lies on the 
inside of the surface, K = -2 tt; and K « 2 tt if P is on the outside of 
the surface. The velocity field may be computed at specified mesh 
points by taking the gradient V p = 

v p - - J J uV[n • V(I)]dS - J J a V(i)dS 
S S 


- J } R,V[n • V(i)]dW + 
W 


(A. 18) 


t 
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To compute the forces and moments due to the body pressure 
distribution, the local velocity at each panel must be determined. 
The following relationships are used: 

V = V* + v (A. 19) 

v = VLZ + VMa + VNn (A. 20) 

where "v" is the perturbation velocity and [Z,m,n] are unit vectors 
in the local panel coordinate system. The normal velocity component 
is obtained directly from the source strength VN * 4 tt a, while the 
tangential components VL and VM are evaluated from the gradient of 
p. The gradient is based on an assumed quadratic doublet 
distribution between three adjacent panels in a given direction 
centered on the panel under consideration. Once the total velocity 
at the center of the panel is known, the force coefficient 
contribution from panel K is 


AF 


K 


, = -c * are^ * V + C * ARE^ * (t k ) 


(A* 21) 


K 


K 


where 


2 2 

C = 1 - v* /V 
P K K . 


(A. 22) 

and the pressure coefficient is evaluated at the panel center. 


is the local skin friction coefficient and is zero during the first 
potential solution; and t^ are unit vectors normal to the* panel's 
mean plane and in the direction of the local velocity 

(t^ = V k/l V kl)’ Now if K 1 and K 2 are the first and last P anels on a 
component (for example, the free-tip), then the component force 

vector is 
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(A. 23) 


F 

-£ = 

q„ 


K 2 

I 

K=K, 


ft 

00 


and the force coefficient vector is computed as 



F /q 
P c 


REF 



AF 

— - /S 
q ' REF 


(A. 24) 


where is a user-specified reference area. The contribution of 

panel K to the moment coefficient vector is 



\ x 

q oo* S REF * L 


(A. 25) 


where is the position vector of the control point on the panel K, 

relative to the moment reference point; for calculation of pitching 
moment coefficient L = CBAR, the reference chord specified by the 
user. Coordinate trajnsf ^are then applied to the force and 
moment coefficient vectors to obtain C^, and C m in wind axes. 
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(4) If M^oCAL ^ ^-*20 everywhere, a check was done to verify 
whether the highest predicted sectional value of the 
tip or inboard section was in the linear range of the 2-D 
experimental lift curve (Figure B.l). Both cases, M^ = 
0.70 and M^ = 0.49, were verified in this manner. 

(5) If nonlinear characteristics were encountered, a*p was 
reduced gradually until maximum values were within the 
linear range for the specified Mach numbers. 

The procedure outlined above is based on the assumption that VSAERO 
predicts the spanwise loading distribution with good accuracy, since 
it is the maximum predicted sectional which is checked against 
the data of Figure B.l. Also, even if some local Mach numbers are 
greater than one, then the calculated sectional lift may still lie 
in the linear range (as long as the supersonic flow region is not 
too large). This is based on the fact that even after sonic flow is 
encountered over an airfoil, M^, may be increased appreciably before 
the onset of nonlinear effects due to compressibility (Reference 
12 ). 

Because of the comparisons with 2-D experimental data, three- 
dimensional viscous effects which may cause separation cannot be 
taken into account using this approach. However, the final 
selection of and a>p are believed to be sufficiently conservative 
to justify its use. 
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Figure B.l: 2-D lift characteristics of the V23010-1.58 airfoil 

in compressible flow (taken from Reference 9) 
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APPENDIX C: 


COMPUTED LIFT AND PITCHING MOMENT COEFFICIENTS 
AT Ho, = 0.70 and ^ - 0.49 

for TIP 2, TIP 3, TIP 4 AND TIP 6 


The data in the following tables follow from the incidence 
angle limits specified in Table 3.3. Observe that generally , the 
predicted growth of the separated flow region agrees with 
corresponding increases in the tip angle of attack. In some cases, 
however, there may be sudden increases in the size of the separated 
flow region, or there may be irregularities in the trends. For 
example, in Table C.l at A0 = -5° (M*, = 0.70), as varies from -1° 
to -5°, the separated flow region on the lower surface increases to 
14.9% and then decreases to 9.0%. These discrepancies are 
associated with the fact that these data are initial estimates only 
and do not represent converged solutions (see Section 3.2.1). 

The asterisks indicate the points which were not used to 
calculate local lift and pitching moment slopes. These did not fall 


within the linear range established by the two lowest values of 
or Lift and pitching moment data were also produced for ay ° 

6° to verify the trends in C L ^ and established by the data 
points falling within the limits of Table 3.3; if within the linear 
range, this fourth point was then used to calculate the slope. At 


= 0.7, VSAER0 output showed some local Mach numbers (on the tip 
surface) exceeding 1.2 for = 6°/a.p = 4°, at = 0.70. Also, for 
tip angles of -5° and -8°, the lift deviates from the linear range 
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due to separation, so those points were also not used in the slope 
computations. The final value of the and slopes were 
obtained by taking the average of the slopes for each A0 at the 
given Mach number. When data points in the nonlinear range were 
eliminated, the slopes were always in good agreement with each 
other. 

Note that at = 0.70, the combination of high Reynolds number 

and large compressibility effects (although subcritical) results in 

chordwise pressure gradients that cause larger regions of separated 

flow than at = 0.49. This normal trend further indicates that 

converged solutions in the nonlinear range could probably be 

obtained with further use of the wake stitching procedure. Observe 

from Table 3.1 and the separation patterns of Tables C.1-C.4 that as 

A of the outboard half of the tip decreases (consider tips 2, 3, 

and 6), the size of the separated flow regions also decrease as 

expected. This is a critical observation, since (C ) T would 

a 

decrease at the onset of stall, thereby illustrating the need for 

compromise between high sweep angles [high negative (C ) T at low 

a 

Mach numbers] and a decrease in response [loss of (C ) T at high M^ 

a 

and Reynolds numbers] associated with these higher sweep angles. 
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Table C.l: Predicted Compressible Lift and Pitching Moment for 

TIP 2 


(1) 

• 

o 

It 

8 

22 

70, 

Re# = 

3.390 x 10 6 




“w 

or 


A9 

( C L^T 

<C "0.25 )t 

% Separated Flow 
Upper Lower 

0 

-8 


-8 

-0.5169 

0.1129* 

0.0% 

10.5% 

2 

-6 


-8 

-0.3549 

0.0704* 

0.0% 

9.24% 

4 

-4 


-8 

-0.1616 

0.0271 

1.2% 

9.0% 

6 

-2 


-8 

0.0486 

-0.0194 

— 

— 

0 

-5 


-5 

-0.3589 

0.0691* 

0.0% 

9.0% 

2 

-3 


-5 

-0.1595 

0.0237 

1.4% 

14.9% 

4 

-1 


-5 

0.0558 

-0.0248 

2.2% 

2.4% 

6 

+1 


-5 

0.3030 

-0.0746* 

— 

2.4% 

0 

-2 


-2 

-0.1602 

0.0213 

1.6% 

16.7% 

2 

0 


-2 

0.0576 

-0.0285 

2.4% 

2.2% 

4 

+2 


-2 

0.3063 

-0.0786* 

6.6% 

1.5% 

6 

+4 


-2 

0.5532 

-0.1295* 

— 

— 

(2) 

© 

II 

X 

491 

, Re# 

= 2.378 x 10 6 




“W 

CXr£ 


A0 

(Cl)t 

( C m>T 

% Separated Flow 
Upper Lower 

6 

-8 


-14 

-0.3556 

0.0662* 


__ 

8 

-6 


-14 

-0.1784 

0.0345* 

0.0% 

13.8% 

10 

-4 


-14 

0.0031 

-0.00048 

0.0% 

11.2% 

12 

-2 


-14 

0.1878 

-0.0370 

1.4% 

0.75% 

6 

-5 


-11 

-0.1604 

0.0267 



8 

-3 


-11 

0.0238 

-0.0098 

0.96% 

13.5% 

10 

-1 


-11 

0.2095 

-0.0471 

1.73% 

0.87% 

12 

+1 


-11 

0.3991 

-0.0856 

2.33% 

0.0% 

6 

-2 


-11 

0.0357 

-0.0158 



— 

8 

0 


-8 

0.2201 

-0.0530 

1.8% 

0.27% 

10 

+2 


-8 

0.4139 

-0.0933 

2.7% 

0.0% 

12 

+4 


-8 

0.6150 

-0.1357 

7.5% 

0.0% 
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Table C.2: Predicted Compressible Lift and Pitching Moment 

for TIP 3 


(1) 

°w 

Mo, = °' 70 » Re# = 

ot*p A 9 

3.390 x 10 6 

(C L ) T 

(C "0.25 >T 

% Separated Flow 
Upper Lower 

0 

-8 

-8 

-0.5175 

0.0974* 

0.0% 

8.5% 

2 

-6 

-8 

-0.3583 

0.0603* 

0.0% 

6.8% 

4 

-4 

-8 

-0.1639 

0.0210 

0.0% 

12.4% 

6 

-2 

-8 

0.0474 

-0.0207 

— 

— 

0 

-5 

-5 

-0.3610 

0.0588* 

0.0% 

6.6% 

2 

-3 

-5 

-0.1611 

0.0183 

0.0% 

13.4% 

4 

-1 

-5 

0.0544 

-0.0248 

1.3% 

1.6% 

6 

+1 

-5 

0.3034 

-0.0686* 

— 

— 

0 

-2 

-2 

-0.1611 

0.0166 

0.0% 

2.6% 

2 

0 

-2 

0.0571 

-0.0277 

1.6% 

1.1% 

4 

+2 

-2 

0.2074 

-0.0716* 

2.5% 

0.0% 

6 

+4 

-2 

0.5576 

-0.1174* 



(2) 

M_ = 0.491, Re# 
00 7 

= 2.378 x 10 6 




“w 

Op 

A0 

(c L ) T 

( C m^T 

% Separated Flow 






Upper 

Lower 

6 

-8 

-14 

-0.3600 

0.0556* 



M Iia 

8 

-6 

-14 

-0.1828 

0.0276 

0.0% 

9.9% 

10 

-4 

-14 

-0.0019 

-0.0026 

0.0% 

3.9% 

12 

-2 

-14 

0.1855 

-0.0343 

0.0% 

0.0% 

6 

-5 

-11 

-0.1643 

0.0206 

— 



8 

-3 

-11 

0.0206 

-0.0108 

0.0% 

0.0% 

10 

-1 

-11 

0.2076 

-0.0436 

0.0% 

0.0% 

12 

+1 

-11 

0.3975 

-0.0767 

1.1% 

0.0% 

6 

-2 

-8 

0.0338 

-0.0164 





8 

0 

-8 

0.2191 

-0.0488 

0.0% 

0.0% 

10 

+2 

-8 

0.4132 

-0.0838 

0.0% 

1.5% 

12 

+4 

-8 

0.6188 

-0.1226* 

2.7% 

0.0% 
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Table C.3: Predicted Compressible Lift and Pitching Moment 

for TIP 4 


(1) 

= 0. 

70, 

Re# = 

3.390 x 10 6 

- 



“w 

Op 


A0 

(C l >t 

<C "0.25 )t 

% Separated Flow 
Upper Lower 

0 

-8 


-8 

-0.5104 

0.0800* 

0.0 

6.21% 

2 

-6 


-8 

-0.3526 

0.04721 

0.0 

3.55% 

4 

-4 


-8 

-0.1630 

0.0142 

0.0 

3.47% 

6 

-2 


-8 

0.0454 

-0.0210 

— 

— 

0 

-5 


-5 

-0.3562 

0.0469* 

0.0 

2.96% 

2 

-3 


-5 

-0.1604 

0.0125 

0.0 

1.89% 

4 

-1 


-5 

0.0536 

-0.0251 

0.0 

0.0 

6 

+1 


-5 

0.3043 

-0.0633* 

— 

— 

0 

-2 


-2 

-0.1601 

0.0116 

0.0 

0.91% 

2 

0 


-2 

0.0564 

-0.0269 

0.0 

0.0 

4 

+2 


-2 

0.3093 

-0.0656* 

1.6% 

0.0 

6 

+4 


-2 

0.5595 

-0.1052* 


— 

(2) 

o 

1! 

s 8 

491 

, Re# 

= 2.378 x 10 6 

- 



“w 

Ct»p 


A0 

( C L^T 

( C m)T 

% Separated Flow 
Upper Lower 

6 

-8 


-14 

-0.3556 

0.0419* 

_ 

__ 

8 

-6 


-14 

-0.18452 

0.0200* 

0.0 

12.1% 

10 

-4 


-14 

-0.0008 

-0.0058* 

0.0 

0.0 

12 

-2 


-14 

0.1755 

-0.0300* 

0.0 

0.0 

6 

-5 


-11 

-0.1675 

0.0145* 





8 

-3 


-11 

0.0132 

-0.0108 

0.0 

0.0 

10 

-1 


-11 

0.1993 

-0.0385 

0.0 

0.0 

12 

+1 


-11 

0.3881 

-0.0669 

0.0 

0.0 

6 

-2 


-8 

0.0324 

-0.0173 

— 

__ 

8 

0 


-8 

0.2189 

-0.0455 

0.0 

0.0 

10 

+2 


-8 

0.4127 

-0.0758 

0.0 

0.0 

12 

+4 


-8 

0.6198 

-0.1095 

0.99% 

0.0 
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Table C.4: Predicted Compressible Lift and Pitching Moment 

for TIP 6 


(1) M*, = 0.70, Re# = 3.390 x 10 6 


°w 

“T 

A0 

(C l >t 

(Cm 0.25 )T 

% Separated Flow 
Upper Lower 

0 

-8 

-8 

-0.5069 

0.0715* 

0.0% 

2.0% 

2 

-6 

-8 

-0.3514 

0.0409 

0.0% 

1.0% 

4 

-4 

-8 

-0.1614 

0.0096 

0.0% 

0.0% 

6 

-2 

-8 

0.0445 

-0.0225 

— 

— 

0 

-5 

-5 

-0.3537 

0.0407* 

0.0% 

0.8% 

2 

-3 

-5 

-0.1586 

0.0085 

0.0% 

0.0% 

4 

-1 

-5 

0.0516 

-0.0254 

0.0% 

0.0% 

6 

+1 

-5 

0.2986 

-0.0592* 

— 

— 

0 

-2 

-2 

-0.1575 

0.0080 

0.0% 

0.0% 

2 

0 

-2 

0.0546 

-0.0267 

0.0% 

0.0% 

4 

+2 

-2 

0.3050 

-0.0612 

0.0% 

0.0% 

6 

+4 

-2 

0.5521 

-0.0969* 



(2) 

o 

II 

8 

491, Re# = 

2.378 x 10 6 




°W 

“T 

A0 

(C L )x 

( C m>T 

% Separated Flow 
Upper Lower 

6 

-8 

-14 

-0.3561 

0.0352* 

__ 

- 

8 

-6 

-14 

-0.1867 

0.0153* 

0.0% 

0.0% 

10 

-4 

-14 

-0.0057 

-0.0093* 

0.0% 

0.0% 

12 

-2 

-14 

0.1726 

-0.0296* 

0.0% 

0.0% 

6 

-5 

-11 

-0.1666 

0.0095* 

— 


8 

-3 

-11 

0.0137 

-0.0132 

0.0% 

0.0% 

10 

-1 

-11 

0.1966 

-0.0375 

0.0% 

0.0% 

12 

+1 

-11 

0.3800 

-0.0617 

0.0% 

0.0% 

6 

-2 

-8 

0.0296 

-0.0179 




8 

0 

-8 

0.2189 

-0.0444 

0.0% 

0.0% 

10 

+2 

-8 

0.4075 

-0.0706 

0.0% 

0.0% 

12 

+4 

-8 

0.6076 

-0.0995* 

0.0% 

0.0% 
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